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Recently the problem of the maximum r ss limit of a stable neutron 
star has drawn the attention of many people, because it may seriously 
affect the models of pulsars and other phenomena which are most likely caused 
by the presence of neutron stars. In this brief report I wish to compare 
the neutron star models by various groups, including the most recent results 
I am aware of, those by the Cornell group (Boozer-Salpeter ) and -chose by 
the Kyoto group (s • Ikeuchi, T. Mizutani, S. Nagata, R. Tamagaki, >nd 
C, Hay ash JJ, 

In Figure 1, the curve marked IDEAL represents the models with no 
nuclear interactions, originally constructed by Oppenheimer and Volkoff. 1 
The models by Harrison, Thorne, Wakano and Wheeler 2 approxin itely lie on 
the same curve in the neutron star region. The dotted curves marked 
and are the models constructed by Tsuruta and Cameron 3 and subsequently 
used by Thorne and others at CIT 4 in their calculations of the moments of 


inertia of neutron stars, etc. I decided to show these curves also, because 


these results (especially the models j have been used by Ostriker and 
others in their pulsar studies. The solid curves show neutron star models 
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constructed after the pulsar discovery. The curve (.1) is by Cohen et al. 8 , 
(2) is by the Cornell group (Booze r-Salpe ter ) , (3> is by Wang et al. 7 , 
and the curves (4) through (J) are obtained by the Kyoto group. 

In all models shown here except the first (those marked IDEAL), nuclear 

interactions are taken into account. In the models marked as V Q and V . 

3 Y 

the nuclear potentials of V Q and V type, respectively, by Lovinger and 
Simmons 8 are used. In the models (l), the modified Lev inger- Simmons neutron 
gas models are used. Boozer and Salpeter (the curve (2)) made use of the 
neutron gas calculations (2b) by Nemeth and Sprung^ with the soft-core Reid 
potentials. Wang et al. (the curve (3)) took the average of the soft-core 
Reid potentials and several other potentials, but they all give the similar 
equations of state. The results by the Kyoto group are obtained in the 
following way. For densities p less than p Q =- 8 x 10 14 gm/cm 3 where the 
neutron matter calculations from the concept of "nuclear potential" 
becomes unreliable, the one-boson exchange hard-core potentials constructed 
by Kishi were used in the models (4), called OBEP-K, and the one-pion 
exchange potentials with a Gaussian type soft core constructed by Tamagaki 
were used in the models ( 5 ) called OPEG-T. For p > p Q in these models, 
the equation of state obtained in this manner was parabclically extrapolated 
in the logarithmic scale. In the models (6) and (j), the equation of 
state P = e/3 and P = e, respectively, (where e is the energy density) 
were used for p > =* 5 x 10 13 gm/cm 3 where the Zel' dovich type scalar 

or vector interactions are assumed to become applicable, the method (4) 
was used for p < p , and the intermediate regions have been interpolated. 
Similar results are obtained if the OPEG-T or the Hamada -Johns ton potentials 
are used for p < p • 
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We see, first of all, that drastically different masses are obtained 

depending on how the inter- particle interaction problem is treated. In 

the models ( 2 ) and ( 3 ) the calculations are terminated where the nuclear 

potential approach is thought to become unreliable. If the calculations 

are extended to higher densities, the curve ( 2 ) reaches the mass peak of 

about 1 M 0 , similar to the curves (4-) and (5) by the Kyoto group. Near 

the breakdown point, p , the models (3), (4-), ( 3 ), ( 6 ) and V have small mass 

u 0 

values of only around 0.2- 0.4- M 0 , while the other models shown in the 
figure have larger masses. At the mass peak, the curve ( 7 ) reaches almost 
3M ® . It may be pointed out that these models obtained by different methods 
nevertheless possess a few common points as mentioned below. As the mass 
increases from the minimum value to approximately 0.2M © , the corresponding 
radius decreases quickly from about 300 km to about .10 km or so. This 
occurs a little below or around the nuclear density, depending on the models. 
For densities of around 10 lb gm/cm 3 the radius is around 10 km, and for 
densities of around 10 1G gm/cm 3 the radius is around 5 Km. Our studies also 
show that the effects of the presence of hyperons and protons on the mass 
limit are much smaller than the effects of nuclear forces, if the interactions 
among all baryons are basically similar. It will be hard to draw definite 
conclusions from the above results, but the following comments may be valid. 

In the regions where the concept of "nuclear potential" is still valid 
(for approximately p < 10 13 gm/cm 3 ), using realistic potentials alone is 
not enough. For instance, quite different results are obtained in the 
models (2; and ( 3 ), even though the both groups used the similar potentials 
(the soft-core Reid potentials). In this respect (of applying realistic 
potentials in a realistic way), I feel that the models ( 3 ), the 0FEG-T, 
are the best recommended (among the models shown here). However, it may 


be noted that many-body interactions are generally neglected and that the 
nuclear potential term is treated non-relativisticaUy in the work done so 
far. The net effect of the corrections to these approximations seems to 
lower the densities and increase the masses near the mass peak. Thus the 
models used already in various pulsar studies seem to be not far from 
reality. 

The exact behavior in the regions above about 10 1Lj gm/cm 3 seems to 
be beyond the knowledge of present-day physics. Since the maximum mass of 
some of the models lie in these high density regions, it is hard to answer 
the question of how high masses stable neutron stars can have. It will 
be fortunate if such a question can be answered rather from the observational 
side. Will it be an impossible dream, if one contemplates that more detailed 
pulsar observations might help the break-through of the difficulties facing 
us today in particle physics and some other fields? 

In the pre-pulsar discovery periods, I did not see much sense in going 
beyond our and models of neutron stars. But, today, with more 
observations of pulsars becoming available, better theoretical work with 
the cooperation of experts in various fields, including the effort to 
construct more realistic neutron star models, seems very desirable. 

I wish to thank Professors Hayashi, Nagata, Salpeter and Tamagaki, and 
Messrs. Boozer, Ikeuchi and Mizutani, for giving me the results of their 
calculations 10 before publications. 
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